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Abstracts 

The two dimensional CP 1 model with 9 term is simulated. We compute the topological 
charge distribution P(Q) by employing the "set method" and "trial function method", 
which are effective in the calculations for very wide range of Q and large volume. The 
distribution P(Q) shows the Gaussian behavior in the small f3 (inverse coupling constant) 
region and deviates from it in the large f3 region. The free energy and its moment are 
calculated as a function of 9. For small (3, the partition function is given by the elliptic 
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theta function, and the distribution of its zeros on the complex 6 plane leads to the first 
order phase transition at 6 = it. In the large f3 region, on the other hand, this first order 
phase transition disappears, but definite conclusion concerning the transition is not reached 
due to large errors. 



1. Introduction 

The two dimensional CP N ~ X model is a suitable laboratory to study dynamics of 
QCD. Topology is expected to play an important role in non-perturbative nature of the 
dynamics of such theories. Recently numerical study of the topological aspects of the 
QpN-i moc j e i mac le much progress [0 @. However a full understanding of the dy- 
namics of the model requires study of an additional contribution of the imaginary part of 
the action, i.e., 9 term. The degeneracy of the different topological sector is resolved into 
a unique vacuum labeled by a parameter 9. As shown by some analytic studies, various 
models with the 9 term, in general, exhibit a rich phase structure H Q |5| 0. It is then 
worthwhile to study effects of the 9 term to the CP N ~ X model 0. From a realistic point 
of view also, it is significant to clarify the matter of strong CP violation in QCD. 

Introduction of the 9 term does not allow ordinary simulations because of the complex 
Boltzmann factor. An idea to circumvent the problem is to introduce the constrained 
updating of the fields, in which the topological charge, being a functional of the dynamical 
fields, is constrained to take a given value Q. So the phase factor e 1 ®® is factored out so 
that the partition function is given by the summation of the probability distribution P(Q) 
weighted by e l9 ® over all possible values of the topological charge Q. This algorithm was 
adopted in simulating the two dimensional U(l) gauge model 0. 

So far topological aspects of the CP 1 model has been studied considerably well both 
theoretically and numerically. Most works are, however, limited to the theory without the 
9 term. This is one of our motives for studying effects of 9 term on the model by means 
of Monte Carlo simulations. We present here the results of P(Q) and the free energy F(9) 
and its moments as a function of 9 by surveying comparatively wide range of Q. Herein the 
two techniques are involved; one is to take the set method [§, and the other is to update by 
modifying the action to the effective action by introducing trial probability distributions in 
each set. These enable one to reach very large Q's. In ref. |]], the model was investigated, 
and the nature of the dilute gas approximation was clarified. In the present paper, we are 
interested in the phase structure in 9 and (3 ( inverse coupling constant ) space. We do 
simulations extensively for various j3 in larger volume V and wider range of Q using the 
techniques. 

From a viewpoint of condensed matter physics as well, the phase structure of the CP 1 
model is worthwhile to study. The antiferromagnetic quantum Heisenberg chain with spin 
s is, in the large spin limit, mapped to the two dimensional 0(3) non-linear sigma model 
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(CP 1 model), as an effective theory describing the low energy dynamics. The topological 



nature appears through the 9 term with 9 = 2irs |T(J, and then its effect is expected to 
distinguish the dynamics between the s = integer and half-odd integer cases. This feature 
is stated as the Haldane conjecture that the s = integer antiferromagnetic Heisenberg chain 



develops gap, while s = half-odd integer one is gapless [11 



The above mapping is based on spin wave approximation in the large s limit. However 
the approximation is believed to be good even for very small s. It is well known that s = 1/2 
antiferromagnetic Heisenberg model is gapless, and correspondingly the non-linear sigma 
model with 9 = n would be critical. There have been analytic arguments assuring this 
||12j| . So far, however, only a few works have been done in terms of numerical calculations 
13[] . The present paper also concerns this issue numerically. 

As shown in this paper, P(Q) shows qualitatively quite different behavior for small 
and large coupling constants; it shows clearly the Gaussian behavior in the small f3 (inverse 
coupling constant) region. We will discuss about the possible first order phase transition 
deducing from the Gaussian distribution. It is based upon analytical discussion by using 
the third elliptic theta function and the Poisson sum formula. In the large (3 region, P(Q) 
systematically deviates from the Gaussian. We show the structural difference of F{&) and 
its moments from the small f3 region at finite 9. Near 9 = n, however, we are not able 
to draw a definite conclusion about the phase structure due to the large errors. We will 
discuss this matter in detail. 

In the following section we fix the notations and present a brief account of the al- 
gorithm of the simulations. In section 3 we give the results. In section 4, the partition 
function zeros are discussed. Conclusions and discussion are presented in section 5. 



2. CP model in two dimensions 

2.1. notations 

We consider the CP 1 model with 9 term on a two space-time dimensional euclidean 
lattice defined by the action 

S e =S - iOQ, 

S = f3 ^ ^ | Zn+^Zn | , ^ ^ 
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where Q is a topological charge, and z a ^ n is a two component complex scalar field (a = 1, 2) 
at site n constrained by 

a 

and couples with the one (z is complex conjugation of z) at the nearest neighbor 

sites n + n(fx = 1, 2). 

The partition function as a function of the coupling constant (3 and 9 is defined by 
Z{9) = I Y[dz n dz n e- Se I I Y[dz n dz n e- S (2.2) 

n n 

and the free energy density F(9) is given by 

F(0) = -ilogZ(0), (2.3) 

where V is volume of the system. The topological charge Q is the number of times the 
fields cover the sphere S 2 . The lattice counterpart we adopt is that of the geometrical 
definition in ref. [Q; the charge density Q(n*) at dual site n* is given by 

Q(n*) = — Im{ In \TrP(n)P(n + l)(P(n + 1 + 2)1 

+ ln[TrP(n)P(n + l + 2)P(n + 2)]|, (2.4) 

where P(n) a p = z an Zf3 n and n is the left corner of the plaquette with center n*. This 
amounts, in terms of z, to the topological charge 

Q = ~. / . tuuifin.u + 9 n -\-u v 9 n -\- U u 9 n v ) 1 (2-5) 

47T z — ' 

where = arg{z n z n+M }. 

In order to simulate the model with the complex Boltzmann factor, we follow the 
Wiese's idea 0. It, in principle, introduces the constrained updating of the fields, in 
which the topological charge, being a functional of the dynamical fields, is constrained 
to take a given value Q. So the phase factor e ld ® is factored out, so that the partition 
function is given by the summation of the probability distribution P(Q) weighted by e 10 ® 
in each Q sector. This amounts, in practice, to calculate first the probability distribution 
P(Q) at 9 = and to be followed by taking the Fourier transform of P(Q) to get the 
partition function Z{9) as 

Z(9)=J2 P (Q)e i6Q , (2-6) 
Q 
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where P(Q) is 

rrr dz dz (<2) e ~ s 

J [[ n dz n dz n e b 

Here JJ n dz n dz n ^' = Yl n dz n dz n 8Q g, i.e. , the integration measure restricted to the 
configurations with given Q, where Sq q is the Kronecker's delta. Note that J^g P(Q) = 1. 
Expectation value of an observable O is given in terms of P(Q) as 

where (O)q is the expectation value of O at 9 = for a given Q sector 

<0>0 = ;n„^*„««e- ' (2 ' 9) 



2.2. algorithm 

We measure the topological charge distribution P(Q) by Monte Carlo simulation by 
the Boltzmann weight exp(— S), where S is defined by (|2.1|) . The standard Metropolis 
method is used to update configurations. To calculate P(Q) effectively, we apply (i) the 
set method and (ii) the trial distribution method simultaneously. In the following, we 
explain briefly the algorithm to make the paper self-contained. All we have to calculate 
P(Q) is to count how many times the configuration of Q is visited by the histogram method. 
The distribution P(Q) could damp very rapidly as \Q\ becomes large. We need to calculate 
the P(Q) at large |Q|'s which would contribute to F{0), (Q)o and (Q 2 )g because they are 
obtained by the Fourier transformation of P(Q) and its derivatives. Further, the error of 
P(Q) at large \Q\ must be suppressed as small as possible. These are reasons why we apply 
two techniques mentioned above. Since P(Q) is analytically shown to be even function 
and is certified by simulation, we restrict to the range of Q to > 0. 

The range of Q is grouped into sets Si ; Si(Q = ~ 3), S^Q = 3 ~ 6), ■•• , 
Si(Q = 3(i — 1) ~ 3z), • • ■ (set method). Monte Carlo updatings are done as follows by 
starting from a configuration within a fixed set Si- When Q of a trial configuration Ct 
stays in one of the bins within Si, the configuration Ct is accepted, and the count of the 
corresponding Q value is increased by one, while when Ct goes out of the set Si, Ct is 
rejected, and the count of Q value of the old configuration is increased by one. This is 
done for all sets Si ; i = 1, 2, • • -. 
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Another of the two techniques is to modify the Boltzmann weight by introducing 
trial distributions Pt(Q) for each set (trial distribution method). This is to remedy P(Q) 
which falls too rapidly even within a set in some cases. We make the counts at Q = 
3(i — l),3(z — 1) + l,3(z — 1) + 2 and 3i in each set Si almost the same. As the trial 
distributions Pt(Q)'s, we apply the form 

Pt(Q) = A t exp[-(a(P)/V)Q% 

where the value of Ci(j3) depends on the set Si, and is a constant. That is, the action 
during updatings is modified to the effective one such as S e g = S + log P t (Q). We adjust 
Ci(0) from short runs to get almost flat distribution at every Q in Si. 

To obtain the normalized distribution P(Q) in the whole range of Q from the counts 
at each set, we make matchings as follows: 

i) . At each set Si ( i = 1, 2, • • • ), the number of counts is multiplied by Pt(Q) at 

each Q. We call the multiplied value Ni(Q), which is hopefully proportional to 
the desired topological charge distribution P(Q). 

ii) . In order to match the values in two neighboring sets S{ and Si+i, we rescale 

N i+1 (Q) so that N l+1 (Q) -> N i+1 (Q) x r, where r =N { (Q = 3i)/N l+1 (Q = 3i), 
the ratio of the number of counts at the right edge of Si to that at the left edge 
of Si+i. These manipulations are performed over all the sets. 

iii) . The rescaled iV$(Q)'s are normalized to obtain P(Q) such that 

P(Q) - W <«> 



3. Numerical Results 

We use square lattices with the periodic boundary conditions. Lattice sizes are V = 
L x L, and L ranges from L = 24, 36, 48 to 72. The total number of counts in each set is 
10 4 . The error analysis is discussed in Appendix. To check the algorithm, we calculated 
the internal energy. It agrees with the analytical results of the strong and weak coupling 
expansions [[I4"| . Using the calculated P(Q), we will estimate the free energy F(6) and its 
derivative (Q)e, respectively. 
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3.1. topological charge distribution P(Q) 

In this subsection we discuss the topological charge distribution P(Q). Partition 
function can be given by the measured P(Q) as in ( |2.6| ) in principle, but we should be 
careful for estimating Z{9) from P(Q). Since P(Q) is very sharply decreasing function 
of Q, its Fourier series Z(9) is drastically affected by statistical fluctuations of P(Q). 
For example, consider two different Q values, say, Qi and Q2 (Qi <C Q%)- Small error 
5P(Qi) at Qi could cause very large effects to Z{9) because P{Qi) itself is sometimes 
much smaller than SP(Qi). So the effort to obtain P(Q) at large Q may be useless if 
we allow these fluctuations at the small value of Q. In order to avoid this problem, we 
first fit the measured P(Q) by the appropriate functions Pfit(Q) an d obtain Z{9) using 
Fourier transforming from Pfit(Q)- We apply the chi-square- fitting to the logarithm of the 
measured P(Q) in the form of polynomial functions of Q 



In the following, we present the results of (3 and volume dependence of P(Q) ■ 

In Fig.l, we show the measured P(Q) for various /3's ((3 = 0.0, 0.5, • • •, 3.5) for a 
fixed volume (L = 24). As (3 varies, P(Q) smoothly changes from strong to weak coupling 
regions. In the strong coupling regions (f3 < 2.0), P(Q) shows Gaussian behavior. In the 
weak coupling regions (2.75 < /3), P(Q) deviates gradually from the Gaussian form, being 
enhanced at large Q compared to the Gaussian. In order to investigate the difference 
between the two regions in detail, we use the chi-square- fitting to logP(Q). Table I shows 
the results of the fittings, i.e., the coefficients a n of the used polynomial J2 n a nQ n for 
various /3's with the resulting x 2 /d.o.f . (i) For (3 < 2.0, P(Q)'s are indeed fitted well by 
the Gaussian form, (ii) For (3 > 2.75, terms up to quartic one are needed for sufficiently 
good fitting. The linear term, in particular, is important for fitting the data at very small Q 
values. The value QMax 5 which is the largest Q of the range in consideration, is also shown 
in the table. It is chosen so that the ratio P(QMax)/P(0) ~ 10 -20 in the weak couplings, 
(iii) Between the strong and weak couplings (2.0 < (3 < 2.75) the fittings according to the 
polynomial turn out to be very poor ( x 2 /d.o.f . ~ 250 ). It may indicate the existence of 
a transitive region between the Gaussian and non-Gaussian regions, (iv) Apart from this 
region, each of the coefficients change smoothly from the strong to weak coupling regions 
as shown in Table I. 




n 
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Here we discuss the volume dependence. In the strong coupling regions, P(Q) is fitted 
very well by Gaussian for all values of V 

P{Q) oc exp (-k v ((3)Q 2 ) . 

where the coefficient «y(/3)(= 02) depends on (3 and V. Fig. 2 shows logKv(P) vs. log V 
for a fixed j3(= 0.5). We see that Ky((3) is clearly proportional to 1/V 

k v (0) = C(J3)/V. 

This 1 /^-dependence of the Gaussian behavior determines the phase structure of the strong 
coupling region. This will be discussed in detail numerically in §3.2 and analytically in §4. 

The proportionality constant C depends on (3. As (3 becomes large, C(/3) monotoni- 
cally increases; C((3 = 0.0) = 10.6, C{(3 = 0.5) = 12.3, C((3 = 1.0) = 15.5. 

Fig. 3. shows the volume-dependence of P(Q) for L = 24, 36, 48, and 72 in the weak 
coupling regions ((3 = 3.0). We do not find the 1/V-law as in the strong coupling regions, 
but a clear volume dependence is observed. It causes the different behavior of F(6) from 
that in the strong coupling regions. 

3.2. Free energy and expectation value of topological charge 

Partition function Z(9) as a function of 9 is given by ( |2.6|) from P(Q). The free energy 

is 

F(O) = -±]0gZ(O). (3.1) 
In general, the n-th order of the moment is given by the derivatives of F{6) 

(Q n )o = -H) n ^r^ ( 3 - 2 ) 

In the strong coupling region, we have seen the Gaussian behavior of P(Q), and the 
1/V-law appears to hold up to L = 72. It is natural to expect that this behavior persists to 
V — ► 00. Let us look at how the 1/V-law affects F{9) and (Q)q. By putting C(0) = 12.3 in 
P(Q) oc exp [- (C((3)/V) Q 2 ] for (3 = 0.5, we calculate F(9) and (Q) from Q and 
Fig's 4 and 5 show their volume dependence. As V is increased, F(9) very rapidly (already 
at L — 6) approaches the quadratic form in 9 from below. Its first moment (Q)e develops 
a peak near 9 = 7r, and the position of the peak quickly approaches tt as V increases. The 
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jump in (Q)e would arise at 9 = tv as V — > oo. It indicates the first order phase transition 
at 9 = tv. 

In the weak coupling regions, on the other hand, we see the different behavior. Fig. 6 
shows F(9) at (3 = 3.0. For 9 < tt/2, F(9) is volume independent, while for 9 > tt/2, the 
clear volume dependence appears, where F(6) decreases as V — > large unlike in the strong 
coupling case. We have checked that the result for L = 20 agrees with that in ref. within 
errors. The expectation value (Q)e is shown in Fig. 7. The singular behavior at 9 = tv, 
which was seen in the small f3 region, disappears. The peak gets round and its locus moves 
towards small 9 as V increases, which is opposite to Fig. 5. 

We should make a remark about the errors in the figures. As a general tendency, 
larger errors arise for larger volume and/or for 9 ~ tv. It is associated with the algorithm 
to calculate Z(6), ( |2.6|) , in which e %e< ^ ~ (— 1)^ for 9 = tv yields large cancellation for 
slowly falling P(Q) (the behavior at large V) in the summation. It causes large errors of 
the observables due to the denominator in (|2.8|) . This is just the same as the so called 
sign problem which is notorious in the quantum Monte Carlo simulations applied to 
systems of strongly correlated electrons. 



4. Gaussian distribution and the partition function zeros 

In the previous sections, we have seen that P(Q) is Gaussian in small f3 region. In this 
section we shall look into the detail of its consequence by paying attention to the partition 
function zeros in the complex £ plane (£ = e %e ). Study of the partition function zeros is 
regarded as an alternative to investigate the critical phenomena. The zeros accumulate 
in infinite volume limit to the critical point, and how fast they approach the point as V 



increases tells the order of the phase transition |16[] [I7|. If Gaussian behavior P(Q) oc 
exp[— (C((3)/V) Q 2 } persists to infinite volume limit, the partition function is expressed by 
the third elliptic theta function 

oo 
Q= — oo 

as 

Z{9)cx# 3 (v, r), 
9 



where p = exp[— C(0)/V] = exp(z7rr) and £ = e %e = exp(z27rz^). In order to look for the 
partition function zeros in the complex £ plane, it is convenient to use infinite product 
expansion of $3 

00 00 

= II (! - P 2 ™) II K 1 + P 2n_1 0(l + V 2n -\- X )\. (4.2) 

m=l n=l 

Zeros of Z(9) are all found easily on the negative real axis of the complex ( plane as 

C = -e~( 2n -V c/v , _ e +(2n-i)c/y 

for n = 1, 2, • • -, 00. In the complex 9 plane, equivalently, these zeros are located at 

9 = n ±i(2n-l)C/V. 

It thus follows that the 1/V-law approaching the critical point 9 C = tt indicates the first 
order phase transition [^] |17|| . 

An alternative to the above way of looking is to use the Poisson sum formula to the 
sum (plf). 

00 00 
Z{9) oc J2 e-°Q 2 / v e ie Q = y^/C ^ e< e ~^ v ^. (4.4) 

Q=—oc n= — 00 

For V ^> 1 and near 9 = n, the sum on the right is well approximated by two terms ( 
n = and 1 ), 

(4.4) » y^/C \e- 02y ^ c + e-(*- 2 -) 2 ^] . (4.5) 

It follows that the partition function has infinite zeros at 9 = n + i(2n + 1)C/V, where 
n is integer. Again the 1/V-law means the existence of the first order phase transition. 
This result is in complete agreement with that from $3 function discussed above. To see 
to what extent the approximation ( f4.5| ) is good, we compare the resulting F{9) and (Q)e 
from (fO|) with those of Monte Carlo simulations. They agree each other. 



5. Conclusions and discussion 

We have seen that P{Q) is Gaussian in the small f3 region. As shown in the last section, 
it leads to the first order phase transition. This behavior is very much like the d = 2 U(l) 
gauge model with 9 term 0, where P(Q) is Gaussian for all values of the coupling constant 
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Tgfl . There the analytic form of P(Q) is given. It may also be interesting to study the 
CP 1 model from the renormalization group point of view, which might show the singular 
behaviors of the renormalization group flows similar to the U(l) case (|]. 

In large (3 region, on the other hand, P(Q) differs from the Gaussian behavior. Con- 
sequently, the free energy F(6) and the moment (Q)e show the quite different behaviors 
from those in the small (3 regions. The signal of the first order transition disappears. To 
understand those behaviors, It would be helpful to consider the dilute gas approximation, 
where instantons of charge Q = ±1 are randomly distributed. Let us assume that the prob- 
ability distribution P n ( P^ ) , in which n instantons (n anti-instantons) generate, obeys the 
Poisson distribution P n = X n e~ x /nl(Pn = X n e~ x /nl). The topological charge distribution 
function P(Q) is given by the modified Bessel's function as P(Q) = e~ x Iq(X), where A 
is average number of instantons (anti-instantons). For A ^> 1, Iq(X) is approximated by 
exp(— Q 2 /2X). The A can then be identified as V/2C, which is natural since the average 
number is proportional to the volume V. As (3 increases, C((3) increases (section 3), that 
is, the average number of instantons decreases; as [3 — > oo (zero temperature limit), the 
configurations vary slowly so that the configurations with large Q are unlikely to contribute 
to the partition function. In large (3 region, the behavior of Iq(X) as a function of Q is 
qualitatively the same with the result of the simulations. Precisely speaking, however, 
they are different, and actually the difference is attributed to the asymptotic scaling of the 
topological susceptibility in ref. |§. 

It is expected from the Haldane conjecture that the second order phase transition 
would occur at 9 = n. We, however, seem to fail confirming it. The first order phase 
transition in small (3 region would have to mutate to the second order one at some /3, if it 
occurred. In the large (3 region, as discussed in section 3, the volume dependence of the 
results is large in the interesting region of 9, and the statistical errors mask the nature. For 
L = 72, the maximal lattice extension of our study, F(9) still changes considerably and 
gets very large errors for 9 > tt/2. Consequently, so do its moments for a wider range of 9. 
This is due to the large correlation length in the large (3 region, and the finite size effect 
is not negligible. The large fluctuations come from the same origin as the so called sign 
problem |[15|| . which arises in the strongly correlated electronic system in the condensed 
matter physics. In order to circumvent the problem, we must address the issue of the 
lattice effect. It is worthwhile to pursue the issue treated in the present paper from the 



the improved point of view such as the perfect action |0|] |2(J . Recently, the second order 
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phase transition has been found numerically by formulating the model in terms of clusters 
with fractional topological charge ±1/2 [ |T5| . 

Some numerical studies of the CP N ~ l model with N > 2 have been done without 
and with the 9 term 0. In the latter case for CP 3 , interestingly, the first order transition 
is observed at finite 9 which is smaller than rr f7]. 
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Appendix 

In this appendix, we discuss briefly the error analysis when "set method" and "trial 
distribution method" are used. 

We consider first the simple case where a single set is adopted and, as trial distribu- 
tion, Pt(Q) = 1- It is known that the counts in the histogram method essentially obeys 
multinomial distribution and that the error of counts at Q (count (Q)) is estimated by the 
variance of the distribution ||. For each Q, the variance is 

a 2(Q) = N . count(Q) / count(Q) 



N \ N 
where N is the total counts. Therefore P(Q) is estimated by 

where SP(Q) = o~(Q)/N. The relative error (5P/P) at large Q is given by 

SP(Q) _ a(Q) = 1 
P(Q) ~ count(Q) y/count(Q)' 



(A.l) 



It could become very large at large Q when P(Q) is rapidly decreasing function of Q. 

When the above two methods are adopted, the relative error decreases as follows. The 
trial distribution method makes count(Q) almost independent of Q. The variance o~(Q) 
also becomes almost constant at each Q. Accordingly, P(Q) is given by 

P(Q) = P t (Q) (count(Q) ± a(Q)) , 

which leads to the relative errors at any Q 

SP(Q) a 



P(Q) count 



constant. 



This is quite an improvement compared to (A.l). When the set method is further used, 
the constant errors do not propagate over different sets H. 
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Table caption 



Table I. The results of chi-square-fitting to log P(Q) in terms of the polynomial J2 n a nQ n for 
various (3. Fittings are performed to the data in the range from Q = to Qiviax- The 
resulting x 2 /do./.'s are also listed. For the data (3 = 0.0, 0.5 and 1.0, Gaussian fitting 
is performed. 



Figure captions 

Figure 1. The topological charge distribution P(Q) vs. Q 2 for (3 = 0.0 to 3.5. The lattice size 
is L = 24. The data only for Q < 21 are plotted. The lines are shown for the guide of 
eyes. 

Figure 2. log<22(= logKy) vs. logV. (3 = 0.5. The 1/V behavior is clearly seen. 
Figure 3. P(Q) vs. Q 2 for (3 = 3.0. The lattice size L is taken to be 24, 36, 48 and 72. 
Figure 4. Free energy F(9) for (3 = 0.5. Lines are shown for V = 16, 25 and 36 in order from 
below. 

Figure 5. The expectation value of the topological charge (Q)e for (3 = 0.5. Lines are shown 
for V = 16, 25 and 36 in order from below. The peak of the curve becomes sharper 
quickly as 9 — > 7r. 

Figure 6. F(6) for (3 = 3.0. L is chosen to be 24 (square), 36 (triangle), and 48 (circle). Values 
of F(9) are plotted based on the parameters a n obtained by the fittings explained in 
the text. The parameters a n for L = 24 are shown in Table I. Those for L = 36 and 
48 are obtained in the same process as for L = 24. The lines are shown for the guide 
of eyes. The volume dependence appears clearly at 9 > n/2. 

Figure 7. (Q)q for (3 = 3.0. L is the same as those in Fig. 6. Error bars for the data of L = 48 
are not drawn because they are too large. 
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Tabic I 



p 


Qlvlax 


x'/d.o.f. 


a 


Oi 


0-2 


13 


a 4 


0.0 


21 


0.91 


-2.57(7) 




-0.018(5) 






0.5 


30 


4.36 


-2.52(7) 




-0.021(3) 






1.0 


21 


1.22 


-2.38(7) 




-0.027(5) 






2.0 


21 


1.15 


-1.97(1) 


-0.012(1) 


-0.058(2) 


0.000(1) 


0.00001(3) 


2.75 


21 


3.75 


-1.21(1) 


-0.20(1) 


-0.191(2) 


0.0064(1) 


-0.00010(3) 


3.0 


16 


8.12 


-0.72(2) 


-0.65(1) 


-0.269(4) 


0.013(4) 


-0.00029(1) 


3.25 


13 


4.92 


-0.33(1) 


-1.55(2) 


-0.257(6) 


0.013(7) 


-0.00030(3) 


3.5 


11 


4.94 


-0.12(1) 


-2.55(2) 


-0.30(1) 


0.021(1) 


-0.00066(7) 
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Figure 1 
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Figure 2 
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Figure 6 
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Figure 7 
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